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1. Introduction 



Exclusive decays of B-mesons provide important tools to test the Standard Model and to 
search for physics beyond it. Hadronic inputs encoding soft physics are not only form factors 
but also light-cone distribution amplitudes (LCDAs). In particular the B-meson LCDAs enter the 
parametrization of the hard-scattering part of hadronic matrix elements of bilocal current opera- 
tors where large momentum is tranferred to the soft spectator quark [1-10]. Impressive progress 
has been made in the calculation of the hard scattering amplitudes entering factorization theorems, 
see e.g. [11-15] for the B —>■ PP case, but one limiting factor for the extraction of fundamental 
parameters is the uncertainty coming from the hadronic input. Recent years have seen several anal- 
yses concerning the renormalization properties [16, [l7|, 18] and the shape of the B-meson LCDAs 
fl 13.' 12.' H3> Hi- UP to now these analyses were restricted to the two-particle case or to 



leading order with the exception of [23]. Here we present the results of [24] for the renormalization 
of the two-particle B-meson LCDAs taking into account mixing with three-parton LCDAs and in 
section (2.3) the results of a new calculation for the combination of three-particle LCDAs *Fa — ^v 
entering the equations of motion. 



2. One-loop calculation with three-parton external state 



The relevant two- and three-parton distribution amplitudes are defined as B to vacuum matrix- 
elements of a non-local heavy-to-light operator, which reads in the two-particle case [ffl] 



(0|^(z)[z,0P v )«(0)|B(p)) 



JbQi) 



(H-^|<M0 + |[<M0-<M0])b 



(2.1) 



a/3 



and in the three-particle case [ |21| ] (the most general decomposition without contraction with a 
light-like vector is given in p5|]): 

(0\q p (z) [z, uz]gG llv {uz)z v [uz, 0] (ftv)a(O) \B(p)) (2.2) 
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We use light-like vectors n± so that every vector can be decomposed as 



in 



(n+-q)—^- + (n--q)—^-+q±^=q+- 
7i 2 = n + -ri-= 2 v = {n + +n_)/2. 



(2.3) 



The computation of the renormalisation properties of the distribution amplitudes requires us to 
consider matrix elements of the relevant operators 

1 



H A(o) 



2% 
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dte mt (0\q(z)[z,0]ri+rh v (0)\H) 



H _((0) = ^ I dte m {0\q(z)[z,0}rt-rh v (0)\H) 



(2.4) 
(2.5) 
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— [5 (to — fc + - <7_) 
1^ 



8((0-k + )]vrlJTT a u 




B:-g s —8(a-k+) vri-TT a u C:g S f . 8 (ca - k+ - q+) U(lt + d)H-TT a u 
v-q [k + q) 

Figure 1: The three leading-order contributions to the matrix element of 0± with a three-parton external 
state. 



with z parallel to n + , 
direction: 



i.e. 



= trij 



v-z 



■ Z-/2 and the path-ordered exponential in the n„ 



[z,0] = Pexp 



ig s / dy^iy) 



= 1 + ig s j^ da z^(az) - g 2 s J Q da J" dpz^ z v (oz) A v (j3z) + . . 
The Fourier-transforms of the different distribution amplitudes are then defined via 

dt I dute^ a+u ^'F(t,u) 



(2.7) 
(2.8) 



0±(W) = 2i J dte ' m ^ = 7^ 



(2.9) 



where F = Wv^a^a^a- Since the renormalization of the operators is independent of the infrared 
properties of the matrix-elements, we can choose an on-shell partonic external state consisting of a 



light quark, a heavy quark and a gluon in equation (2.6). The resulting leading-order diagrams are 
shown in figures [j] and ^. Next-to-leading order (NLO) diagrams are obtained by adding a gluon 



or a quark loop (a ghost loop) in all possible places (for a complete list of diagrams, see [24]). 



Since the operators give rise to c) -distributions in the +-component of the momenta, we chose to 
proceed via the theorem of residues. To be more explicit, we decomposed the loop momentum / 
in light-cone components, picked up the poles in the /--integral and performed the /^-integration 
in dimensional regularization with D = 2 — 2e dimensions. Additional ^-poles arise through the 
l + -integration for diagrams where a gluon is exchanged between the Wilson-line from the operator 
and the heavy-quark field. These are related to the cusp anomalous dimension, see e.g ||2^, p7| ], 
stemming from the intersection of one light-like Wilson line from the path ordered exponential in 
the operator and one time-like Wilson line from the interaction of soft gluons with the heavy quark. 
The additional poles give rise to -r-terms as well as Sudakov logarithms. 
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A-3n ■ ig s {q+e^ l -q ll £ + ) vTT a u8((0 -k+)8(£ -q+) 



Figure 2: Leading-order contribution to the matrix element of with a three-parton external state. 



2.1 + -case 



In the + -case there is no mixing from three-particle distribution amplitudes. y + 3 = at order 
a s . We confirm the result for the anomalous-dimension matrix found in JT6[ ] 

riV^AO = (r^piogii + ^ «(« - ai') -r^a, f + ^.io) 

+ V y CO J v \CO'(co'-CO) co(co — co')J, 



f(a, co') = f(co, co') - 8(co - co') I dco'f(co, co') r« p =4 y 



with 



2.2 0_-case 

The 0_ case is more involved. After coupling-constant and external leg renormalization there 
remains a genuine three-particle term. The renormalization group equation to order a s can be 
written as 



4k 



J dco' y^(co,co';ii)^(co';ix) 



+ J dco'd^Y W 3 (co,co' '^';n)[^ A -^v](co'^';n)] (2.11) 



where yi 1 ' is the result from [TT 



and yi 1} 3 from @] 



yi^O), =4 



yi 1 )( W , W ^) = 7i 1) -ra^-^ 



co' 



CO' 



co' + ^'-co v ; (co' + ^O 



(2.12) 



C A 



®{co' + E'-co) 1 , . ^, . 



where we defined the +-distribution with three variables as 

>(<*>, ffl'.Ol =f(co,co'^')-8(co-co'-^) I dcof(co,co'^' 



(2.13) 



(2.14) 
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2.3 *Pa — *1V -renormalization and equation of motion constraints 

Here we report on an up to now unpublished calculation of the renormalization of the three- 



particle LCDAs \Fa — *Fy. We project on the relevant distribution amplitudes in equation (2.3) 
using r = y^tl + tl-Y5 (although a instead of 7^ yields the same result). The calculations go along 
the same lines as in the previous two cases, even though there is only one leading-order structure 
(shown in figure ^ and NLO diagrams must have one gluon attached to the vertex in order not to 
vanish trivially. For convenience the result is splitted into Cp- and Ca -colour structures. 



e(co-co') 



CO — CO 1 



+ 



X 



-Q(CO-CO') 



e{co-co'y 




CO -co' 


+ 



e(g-g') 

co &(co' - co)~ 
co' co'-co 
co &(co' 



+ <5(w + £ - 



ft) 



ft)' co'-co 



a-? 



+ 5(a) + § - - 50 g /(<p/ + g { -^(ffl' + «' - - w ') 

- -^(w' + 2^' - fti)0(ft)' - co)e(co) + %{co- " <»)®(<») 

+ ^-^( w '+^- w )0(ft,-^')0(^) 



(2.15) 



yg A (o),£, = y^{coM^m-^ ! ) + y^(co^,co ! ^ 



J 1 ) 



.(1) 



ft) 



€2e( *- § W£^iz^e(„,) 



ft)' (co + §) 2 



t; co + t; 



J, co-l 



(2.16) 



with yj 1 ' the same as in equation ( 2.10| ) and y^ defined as in (2.11) with obvious changes. Part of 
this calculation, namely the light-quark-gluon part, has been calculated in a different context and a 



i 1 ) 



different scheme, e.g. in J29|, |28p. 

In two equations from the light- and heavy-quark equations of motion were derived 

ft)0i(ft);ju) + 0+(ft);ju) = /(ft);/i), (ft) - 2A)0 + (ft);/i) + &)</>_ (ft); /i) = J(co;pi), (2.17) 

where I(J)(co; /i) are integro-differential expressions involving the three-particle LCDAs — *Fy 
(*Fa +-Xa an d *Py) respectively. While the second equation was shown not to hold beyond leading 



order in [17, 23] we have checked that the first one is valid once renormalization is taken into 



account. Taking the derivative of the first equation with respect to log/x, inserting 



and using the relation from [17] 



d n dco' m , (i) 

dCO Jo Tl 



(2.18) 



(2.19) 
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one arrives at the following equation 

co^Jdco'd^y%{co,co'^) (¥ A Kr;M)-*vK,<r;M)) 
+ 2 J dco'i l \co,co';^ J Q a dp £ ^-^(W A (p,^n)-W v (p,^)) 

= 2jda/d?±£dp jT ^^ ) (p,«,id , 5';m)(Va(id',« , ;m)-^(<d , ) « , ;m)), 

(2.20) 

which can be proven to hold at order a s by simple insertion of the respective evolution kernels 
( 2.10 ), ( 2.13 ), ( 2.15 ), ( 2.16 ). This non-trivial outcome gives us further confidence concerning the 
renormalization group properties of the LCDAs. 

3. Conclusions 

The presence of 8(co — ft)') log (/i/ ft)) in the renormalization matrices gives rise to a radiative 
tail falling off like (logft))/ft) for large ft). Therefore non-negative moments of the LCDAs are not 
well defined and have to be considered with an ultraviolett cut-off. [TJ, 17, 22, 23] 



r^uv 

ft) )±(ju) = / dcoco 0±(ft);ju' 
Jo 



lim f AUV d(0(0 N z { l\((0,(o',^') = AT = 0,1, 



1 



= 27 7 - 



For 0_ it is interesting to examine the limit 

rMiv 
lim / 

A, 

which is relevant for the calculation of the three-particle contributions to the moments: 

fAf/v ci { f rA-w 



(3.1) 



(3.2) 



y <ift)ft) 0_(<b;ju) = 1 + ^ ( y da' <j)-((o') J dcoco z_ {co,co'\fx) 

- J dco'd$'(2 - D) [¥ A - ¥ v ] § ') / ^ do) a^zL 1 ^ (<», a', § 



(3.3) 



Therefore as stated in [17] three-particle distribution amplitudes give only subleading contribution 



to the first two moments (N = 0, 1) and we have explicitly checked that this statement cannot be 
extended to higher moments (N > 2). 

The next step consists in using the renormalization properties as a guide to go beyond the ex- 
isting models derived from a leading-order sum-rule calculation resulting in *P A = *Fy j j20| ] and to 
analyze their influence on 0_ . Finally, for practical calculations involving three-particle contribu- 
tions, one would need the evolution kernels for the all relevant LCDAs, which will be the subject 
of a future work. 
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